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Abstract 

We construct a time-dependent solution in vacuum string field theory and in- 
vestigate whether the solution can be regarded as a rolling tachyon solution. First, 
compactifying one space direction on a circle of radius R, we construct a space- 
dependent solution given as an infinite number of *-products of a string field with 
center-of-mass momentum dependence of the form e'^"^ Our time-dependent 
solution is obtained by an inverse Wick rotation of the compactified space direc- 
tion. We focus on one particular component field of the solution, which takes the 
form of the partition function of a Coulomb system on a circle with temperature 
E?. Analyzing this component field both analytically and numerically using Monte 
Carlo simulation, we find that the parameter h in the solution must be set equal to 
zero for the solution to approach a finite value in the large time limit — > cxd. We 
also explore the possibility that the self-dual radius R = ^/a' is a phase transition 
point of our Coulomb system. 
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1 Introduction 



The rolling tachyon process represents the decay of unstable D-branes in bosonic and su- 
perstring theories mi21- This process is described in the hmit of vanishing string couphng 
constant by an exactly solvable boundary conformal field theory (BCFT). Study of this process 
has recently evolved into various interesting physics including open-closed duality at the tree 
level, a new understanding of c = 1 matrix theory and Liouville field theory, and the rolling 
tachyon cosmology (see [3H1|3] and the references therein). However, there still remain many 
problems left unresolved; in particular, the closed string emission and its back-reaction [H]. 

One may think that such problems can be analyzed using string field theory (SFT), which 
is a candidate of nonperturbative formulation of string theory and has played critical roles in 
the study of static properties of tachyon condensation (see [ZllHI and the references therein). 
However, SFT has not been successfully applied to the time-dependent rolling tachyon process. 
The main reason is that no satisfactory classical solution representing the rolling process 
has been known in SFT, though there have appeared a number of approaches toward the 
construction of the solutions P CDl UH 1121 UHl HH UH] • Among such approaches, refs. jHl E] 
examined time-dependent solutions in cubic string field theory (CSFT) by truncating the 
string field to a few lower mass component fields and expanding them in terms of the modes 

gnxO ^ 0,±1,±2, ■ ■■)• 

Let us summarize the result of our previous paper (we use the unit of a' = 1). We 
expanded the tachyon component field t{x^) as 

oo 

t(x°) = ^t„coshna;°, (1.1) 

n=0 

and solved the equation of motion for the coefficients t„ numerically (by treating ti as a free 
parameter of the solution). Our analysis shows that the n-dependence of tn is given by 

t„ ~ A-"' (-/?)" , (1.2) 

up to a complicated subleading n-dependence. Here, A is a constant 3^/72^ and p IS a 
parameter related to ti. From the effective field theory analysis, the rolling tachyon solution 
is expected to approach the stable non-perturbative vacuum at large time —>■ oo jJT] . If tn 
behaves like p.2|) . however, the profile of the tachyon field t{x^) cannot be such a desirable 
one: it oscillates with rapidly growing amplitude (see also ()3.1|) and ()3.2|) ): 

t(x°) ~ e("")'/(''"^) X (oscillating term) . (1.3) 

Since the radius of convergence with respect to x^ of the series (jl.ip is infinite for t„ of (jl.2p . 
we cannot expect that analytic continuation gives another t{x^) which converges to a constant 
as x° oo. 
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In order for the series to reproduce a desirable profile, it is absolutely necessary that 
the fast dumping factor A~"^ of (jl.2p disappears. If this were the case and, in addition, if tn 
were exactly given by 

tn = {-f3r, (1-4) 

analytic continuation of the series would lead to 

t(x°) = -l + \-^ + 1 0, (1-5) 

which approaches monotonically a constant as x° oo. This particular t(x°) has another 
desirable feature that it becomes independent of x° when j3 = and 1, which may correspond to 
sitting on the unstable vacuum and the stable one, respectively. Since CSFT should reproduce 
the rolling tachyon process, it is expected that the behavior ()1.2|) is an artifact of truncating 
the string field to lower mass component fields and that some kind of more sensible analysis 
would effectively realize A = 1. 

The purpose of this paper is to study the rolling tachyon solution in vacuum string field 
theory ( VSFT) jTHl HH 1201 HI] , which has been proposed as a candidate SFT expanded around 
the stable tachyon vacuum. The action of VSFT is simply given by that of CSFT with the 
BRST operator Qb in the kinetic term replaced by another operator Q consisting only of 
ghost oscillators. Owing to the purely ghost nature of Q, the classical equation of motion 
of VSFT is factorized into the matter part and the ghost one, each of which can be solved 
analytically to give static solutions representing Dp-branes. In fact, analysis of the fiuctuation 
modes around the solution has successfully reproduced the open string spectrum at the unsta- 
ble vacuum although there still remain problems concerning the energy density of the solution 
1221 123 1211 ■'■f ^6 '^^^ similarly construct time- dependent solutions in VSFT without trunca- 
tion of the string field, we could study more reliably whether SFT can reproduce the rolling 
tachyon processes, and furthermore, the unresolved problems mentioned at the beginning of 
this section. 

Our strategy of constructing a time-dependent solution in VSFT is as follows. First we 
prepare a lump solution depending on one space direction which is compactified on a circle 
of radius R. Then, we inverse- Wick-rotate this space direction to obtain a time-dependent 
solution following the BCFT approach pp. The lump solution of VSFT localized in uncom- 
pactified directions has been constructed in the oscillator formalism by introducing the cre- 
ation/annihilation operators for the zero-mode in this direction In the compactified case, 
however, we cannot directly apply this method. We instead construct the matter part $™ of 
a lump solution as an infinite number of *-products of a string field Qj,; $™ = Qh*Qh* ■ ■ ■*Qb 
(the ghost part is the same as that in the static solutions). Since the equation of motion of $™ 
is simply <l>™ * $™ = $™, this gives a solution if the limit of an infinite number of *-product 

"'^See also jSHlEEl for recent attempts to this problem. 
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exists j^Zj. As the constituent fit, we adopt the one which is the oscillator vacuum with respect 
to the non-zero modes and has the Gaussian dependence e"''^^/^ on the zero-mode momentum 
p in the compactified direction. Finally, our time-dependent solution is obtained by making 
the inverse Wick rotation of the compactified direction X —>■ —iX^ or —iX^ + nR. 

After constructing a time-dependent solution in VSFT, our next task is to examine whether 
it represents the rolling tachyon process. Our solution consists of an infinite number of string 
states, and we focus on one particular component field (we adopt the same symbol as 

the tachyon field in the CSFT analysis). This t{x^) has the expansion with coshnx° 

replaced by cosh(r;,x''/-R). What is interesting about t{x^) is that it takes the form of the 
partition function of a statistical system of charges at sites distributed with an equal spacing 
on a unit circle. The temperature of this system is R^. The charges interact through Coulomb 
potential and they also have a self energy depending on the parameter b of Qb- The partition 
function is obtained by summing over the integer value of the charge on each site keeping the 
condition that the total charge be equal to zero. 

What we would like to know about t{x^) are particularly the following two: 

• Whether t{x^) has a profile which converges to a constant as x^ — >■ oo. 

• Whether the critical radius i? = 1 in the BCFT approach |28t l^ is required also in our 
solution. 

That we have to put i? = 1 in our solution is also natural in view of the fact that the correct 
value —1 of the tachyon mass squared is reproduced from the fiuctuation analysis around 
the D25-brane solution of VSFT [221 12ni 121] • For these two problems, we carry out analysis 
using both analytic and numerical methods. In particular, we can apply the Monte Carlo 
simulation since t{x^) is the partition function of a Coulomb system on a circle. We find that 
the coefficient t„ of our VSFT solution has a similar n-dependence to with A depending on 
the parameter b. This implies that the profile of t{x^) is again an unwelcome one for a generic 
value of b: it is an oscillating function of x° with growing amplitude. However, we can realize 
A = 1 by putting 6 = and taking the number of in $™ = Q^* - ■ ■*Qhto infinity by keeping 
this number even. These properties seems to hold for any value of R. In order to see whether 
R = 1 has a special meaning for our solution, we study the various thermodynamic properties 
of the Coulomb system t{x^). First we argue using a naive free energy analysis that there could 
be a phase transition at temperature R'^ = 1. Below i?^ = 1 only the excitations of neutral 
boundstates of charges are allowed, but above R'^ = 1 excitations of isolated charges dominate 
the partition function. We carry out Monte Carlo study of the internal energy and the specific 
heat of the system, but cannot confirm the existence of this phase transition. However, we 
find that the correlation function of the charges show qualitatively different behaviors between 
the large and small R^ regions when 6 = 0, possibly supporting the existence of the phase 
transition. 
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The rest of this paper is organized as follows. In section 2, first briefly reviewing VSFT 
and its classical solutions representing various Dp-branes, we construct time-dependent so- 
lutions following the strategy mentioned above. In section 3, we investigate the profile of 
the component field t{x^) both analytically and numerically. In section 4, we argue that our 
solution with 6 = could give a rolling tachyon solution. In section 5, we study a possible 
phase transition at i?^ = 1 through various thermodynamic properties of the system. The 
final section (section 6) is devoted to a summary and discussions. In the appendix we present 
a proof concerning the minimum energy configuration of the Coulomb system. 



2 Construction of a time-dependent solution in VSFT 

In this section, we shall construct a time-dependent solution in VSFT. As stated in section 1, 
we first construct a lump solution which is localized in one spatial direction compactified on a 
circle of radius R. This solution is given as an infinite number of ^-products of a string field 
Qb] * * ■ ■ ■ * ^b- Our time-dependent solution is obtained by inverse- Wick-rotating the 
spatial direction to the time one. Throughout this paper, we use the convention a' = 1. 



2.1 D]7-brane solutions in VSFT 



In this subsection, we briefly review the construction of lump solutions in VSFT describing 
various Dj9-branes in the uncompactified space JHI- VSFT is a string field theory around the 
non-perturbative vacuum where there are only closed string states. Its action is written as 
follows using the open string field 



2 3 

1 1 

= -- I Q|^)--o(^|i(^ 12(^1 ^3)012- 



(2.1) 



The BRST operator Q of VSFT consists of only ghost operators, and it has no non-trivial 
cohomology. The three-string vertex | V3 ) represents the mid-point interaction of three strings, 
and it factorizes into the direct product of the matter part and the ghost one. More generically, 
the matter part of the A^-string vertex | ) representing the symmetric mid-point interaction 
of A^-strings (A^ = 3, 4, ■ ■ ■ ) is given by |S1 EH 



26 „ 



r26 



Af-1 
r=0 



N-1 



X exp I 



r,s=0 



00 
n,m=l 



^ '^nO"n Ps^ c) ^00 PrPs 



n=l 



0;Pr )r 



r=Q 



(2.2) 



4 



where | 0; ) is Fock vacuum of the r-th string carrying the center-of-mass momentum (the 
index r specifying the strings runs from to — 1). Here we use the same convention as 
[THl 1201 1^ • o-n'^^ are the matter oscillators of non-zero modes normalized so that their 
commutation relations are 

k''^^ al^"^] = v'^^SnmS'^' , (n, m > 1) . (2.3) 
The coefficients V^^ are called the Neumann coefficients. In particular, Vqq is given by 

Tr(r — .Q^l 

-In 
21n| 



yrs 



^ . n(r - s) 
2 sm 



(2.4) 
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Note that V^^ depends on although we do not write it explicitly. 
The action ()2.ip leads to the equation of motion 

= * ^ . (2.5) 

Assuming that the solution is given as a product of the matter part and the ghost one, 
vl> = vl>™ (g) the equation of motion is reduced to 

^m^^m^m^m^ (2.6) 

Q^s = -^g*s^s, (2.7) 

where *™ (*^) is the *-product in the matter (ghost) sector. In this paper, we assume that 
the ghost part \E'^ is common to the various solutions, and focus on the matter part equation 

(ESD. 

Classical solutions of ()2.6|) which represent the various Dp-branes in spacetime are given 
in • Let us review the two ways of constructing classical solutions representing the transla- 
tionally invariant D25-brane. One way is to assume that \E'™ is given in the form of a squeezed 
state, the exponential of an oscillator bilinear acting on the vacuum: 



VI/- ) = AT exp I -ir^^, Yl ^-<^<^ ) I 0' ) ' (2-^ 



oo 



m,n=l 



where A/" is a normalization factor. The equation of motion ()2.6|) is reduced to an algebraic 
equation for the infinite dimensional matrix Smn, which, under a certain commutativity as- 
sumption and using the algebraic relations among the Neumann coefficients [30J, can be 
solved to give Smn in terms of 



rrs 

' mn 



s 



CT, T=^(l + X- V(1 + 3X)(1-X)) , (2.9) 
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with the matrices C and X given by 

C^n = i-irSmn, X = CV'\ (2.10) 

Another way is to construct as the shver state j2ZI- Defining the wedge states as 
|iV)o = |0;0)* |0;0)*---* |0;0) = i(0;0|2(0;0|---;v-i(0;0|r^)oi...Ar_i, (2.11) 



N-l 



they satisfy the following property: 

|A^)* |M) = |Ar + M-l). (2.12) 

Taking the limit A^, M — ^ cxo, we have 

|oo)*|oo) = |oo). (2.13) 

Namely, the state | oo ) (sliver state) is a solution to ()2.6|) . It has been proved that the two 
solutions, ()2.8|) and | oo), are identical with each other 

Lump solutions localized in spatial directions can be constructed in the same way as the 
D25-brane solution. Let us denote the directions transverse to the brane by x". In the 
squeezed state construction jTH], we introduce the annihilation and the creation operators for 
the zero-mode in the transverse directions by 

a" = ^p"-^a;", a"^ = — r + ^a;", (2.14) 
°2^v^'°2^v^' ' 

where h is an arbitrary positive constant. Since the zero- modes Oq satisfy the same commuta- 
tion relation (|2.3|) as the non-zero modes, we define the new Fock vacuum | fif, ) by 

aM"|n,) = 0, (n>0). (2.15) 

The new vacuum | Vti, ) with the normalization ( | ) = 1 is expressed in terms of the 
momentum eigenstates as 

I ) = n [jP"" ^'^'^''^'"'^ I 0' ) • (2-16) 

With these oscillators and the coordinate-dependent vacuum | fif, ) , the transverse part of the 
three-string vertex | V3 ) can be written as 



exp 



00 \ 

-2 E E «l::^"'c>i'^^"Mifi.)oi2, (2.17) 

r,s=0,l,2 m,n=0 / 
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in terms of the new coefficients V^J^, which satisfy the same algebraic relations as V^^. There- 
fore, we can construct the Dp-brane solutions just in the same way as the D25-brane solution: 

I ) = I ) ® I ) 

= e^vi-\ V,u f2 S^na>:^aA I 0; p ) ® exp T - ^ S:^^a^a'^A | fi, ) , (2.18) 

^ m,n=l ^ ^ m,n=Q ' 

where the indices /i, v run the directions tangental to the branes (//, z/ = 0, 1, ■ ■ ■ ,25 — p), and 
S'.^^ is given by ()2.9j) with V^^ replaced by V'^^ . This lump solution contains one arbitrary 
parameter 6, the physical meaning of which is not known. It has been shown that the ratio of 
the tensions of Dp-brane solutions is independent of h [33] . Later we will argue that we must 
choose 6 = to obtain a time-dependent solution with the desirable rolling profile. 

Finally, note that, since the modified Neumann coefficients V!^^ satisfy the same algebra 
as the original V^^, the transverse part | \E'™ ) of the lump solution ()2.18|1 can be written as a 
sliver state: 

1*^)= lim i(r]5|2(^^6|---Ar-l(^^6|Viv±)oi...iV-l, (2.19) 

where | V/vx ) is the transverse part of the iV-string vertex. 



2.2 Time-dependent solution in VSFT 

Now let us construct a time- dependent solution in VSFT which possibly represents the process 
of rolling tachyon. This consists of the following two steps: 

• Construction of a lump solution of VSFT localized in one space direction which is com- 
pactified on a circle of radius R. 

• Inverse Wick rotation of the compactified space direction to the time one on this lump 
solution to obtain a time-dependent solution in VSFT. 

Since both the solution and the string vertices have factorized forms with respect to the 
spacetime directions, we shall focus only on this transverse direction of the brane in the rest 
of this paper. 

First, we shall construct a lump solution on a circle. The squeezed state construction 
explained in the previous subsection, however, cannot be directly applied to the compactified 
case since the zero- mode creation/annihilation operators of (j2.14j) are ill-defined due to the 
periodicity ~ x° -|- l-nR. Therefore, we shall adopt the sliver state construction of the lump 
solution. Namely, let us consider 

lim I fib) * • ■ ■ * I 1];,) = lim i(fi;, I ■ ■ ■Ar_i(fi;, I VAr)oi...Ar-l , (2.20) 
^ y N—*oo 
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with a suitably chosen \Qb)- If the hmit — > cxd of (j2.2(J|) exits, it gives a solution of VSFT. 
Taking into account that the momentum zero-mode p in the compactified direction takes 
discrete values p = n/R, we adopt as the state \flb) in ()2.20p the following one which is a 
natural compactified version of ()2.16|) : 



nb)= J2 e~('/')("/^)'|0;n/i?) 



(2.21) 



where \0; n/R) is the momentum eigenstate (and the Fock vacuum of the non-zero modes) with 
the normalization {0;m/R \ 0;n/R) = 6n,m- The A^-string vertex | V)^ ) for the compactified 
direction is given by (j2.2j) with the replacements: 

jdp-^^Y.^ Kp)^R^nfi. \Q]p)^^\Q;n/R). (2.22) 

n 

Vtb ) in the x-representation for the center-of-mass dependence is 



Then the state | fi^, ) * • • ■ * 
given by 

oo oo oo 

no=— ooni=— oo rijv-i=— 
, ^ ^ Af-1 oo ^ N-l 

M- oY: ^n^aial -EE ^^n"^ - ^ E 



oo oo oo 

{x\{ \n,)*...*\n,)) = E E ■■■ E ^ 



exp 



iriQX 



R I 



X 



Af-l 
oo 



n,m=l 



=0 71=1 



00 '^r'^s 



r,s=0 



r=l ' 



(2.23) 



which in the limit N ^ oo should give a lump solution on a circle. In this paper we are 
interested only in the time-dependence of the solution and hence ignore the overall constant 
factor multiplying the solution. 

Our construction of a time-dependent solution of VSFT is completed by making the inverse 
Wick rotation X —>■ —iX^, namely, x — > —ix^ and aj^ — > —iaj^, on this lump solution: 



no=— ooni=— oo ?ijv-i=~c» 
oo A''—! oo 



X exp 

where Qrs is defined by 



^ oo N~l oo ^ Af-l s 

^ E ^ZaWm + . E E ^^°o«n| - ^ E ^-^^^0 I ) ' (2-24) 

n,m=l s=0 n=l r,s=0 ' 



-2 In 



2 sin 



vr r — s 



Qrs — + 2 ('^»':0 ^ 1) — 







21n( 












21n( 









b 

2 



(r = s ^ 0) 
(r = s = 0) 



(2.25) 



Taylor expansion of exp (z J2s J2n ^no'^n^s/-^) ()2.24j) gives an expression of | \E'(x°) ) as 



), where | a ) are the static string states of the 
are the corresponding time-dependent com- 



an infinite summation, | ^{x ) ) = \a) (pa{^ 
form ■ ■ ■a'^ exp(| ^ VmnO'ln'^Dl ^^^d 
ponent fields. In this paper, we shall, for simplicity, focus on the component field of the pure 
squeezed state exp(i ^ ^mn^ln^n) I )• Since we have limAr^oo Vmn — ^mn this component 
field is that for the state representing the unstable vacuum. Denoting this component field by 
t(x°), we have 



"0 ' 



where tno is given by 



"no 



exp 



ni,--- ,njv— 1=— oo 
(no+niH hnjv-i=0) 



-—H{nr;no] 



(2.26) 



(2.27) 



with 



Af-l 



H{nr]no) = ^^Q 



(2.28) 



r,s=0 



Note that the coefficient tng can be regarded as the partition function of a statistical system 
with Hamiltonian H and the temperature R^. In this statistical system, we have charges rir 




Figure 1: There are charges rir on the unit circle. Each charge rir takes integer values, 
and the total charge must be equal to zero. The charges ra^ and rig interact via the Coulomb 
potential Qrs- 

on a unit circle at an equal interval (figure C]). The charges take integer values from — oo to 
+00, and they have the self interaction Qrr and the two-dimensional Coulomb interaction Qrs 
(r 7^ s) between each other. In tno, the charge no at r = is fixed, and there is a constraint 
that the total charge be equal to zero. Note that tno positive definite, tno > 0; is even 
under uq —tiq: 

t—no tno ■ (2.29) 
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On the other hand, t(x°) itself is interpreted as the partition function of the statistical system 
in the presence of the external source / R for uq. 

Solving the constraint ^1,^=0 = to eliminate n^r-i, and are rewritten using 
independent variables without constraint: 

oo / ^ N-2 \ 

^"0 = exp - -— ^ QrsUrUs , (2.30) 



ni,--- ,njv— 2=~cxD 



where Q^s is a (A^ — 1) ^ {N — 1) matrix given by 

Qr,s = Qrs-Qr,N-l-QN~l,s + QN~l,N-l, (r, S = 0, 1, ■ ■ ■ , - 2) . (2.31) 

We have checked numerically that the matrix Qrs is a positive definite matrix. 

In addition to the above | \E'(x°) ) obtained by the inverse Wick rotation X —iX'^, we 
have another time-dependent solution via a different inverse Wick rotation, X — > —iX^ + nR. 
This new solution is obtained simply by inserting (—1)"" into ()2.24|) and ()2.26|) . and satisfies 
the hermiticity condition. As we shall explain in the next section, we expect that this new 
solution with represents the rolling process to the stable tachyon vacuum, while the 

orignal solution given by (j2.24j) represents the rolling in the direction where the potential is 
unbounded from below. 



3 Analysis of the component field t{x^) 

In this section, we study the profile of the component field t{x^) given by ()2.26p both analyti- 
cally and numerically. If our VSFT solution ()2.24|) represents the rolling tachyon solution, the 
component field t{x^) as well as the whole | \Ef(x°) ) should approach zero, namely the tachyon 
vacuum, as x° ^ oo. 

Let us mention the expected riQ-dependence of the coefficient tno, ()2.27p and ()2.30|) . neces- 
sary for t{x^) to have a rolling tachyon profile. Suppose that the no-dependence of tno given 

- — ' 2 

where t„Q has a milder no-dependence than the leading factor e~''"o. We expect that 
lim^^Qo tfi/tn +1 is finite and larger than one, namely, that the series (|2.26p with tno replaced 
by tno a finite radius of convergence with respect to A typical example is tno ~ e"'''"'''. 
Such tno actually appeared in the time-dependent solution in CSFT in the level-truncation 
approximation E] . For this tno ' have 

oo 
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If does not depend on hq, is a periodic function of x° with period 2a, and the 

whole t{x^) cannot have a desired profile: it oscillates with blowing up amplitude e*^^"-*^''^'*"-* as 
x° oo. Even if tng has a mild no- dependence such as t„Q ~ e~'''"'°', it seems very unlikely that 
t{x^) approaches a finite value in the limit x° — oo. These properties persist in the alternating 
sign solution (— l)"°)f:„g obtained by another inverse Wick rotation mentioned at the end of 
section 2. Therefore, it is necessary that the leading term e~°'"'o in ()3.H) is missing, namely, 
we must have a = 0. If this is the case, the series ()2.26|) would have a finite radius of 

convergence, and the analytic continuation would give a globally defined such as ()1.5|) . 
Since tno is positive definite, diverges at the radius of convergence and corresponds to 

the rolling in the direction of the unbounded potential. On the other hand, another t{x^) with 
alternating sign coefficients (—1)"" t„g is expected to be finite at the radius of convergence and 
represent the rolling to the tachyon vacuum. 

In the rest of this section we shall study whether the condition a = is satisfied for the 
present solution. We shall omit the indices of uq and unless confusion occurs. 



3.1 Analysis for i?^ > 1 and i?^ <C 1 

In this subsection, we shall consider the n-dependence of the coefficient t„ of ()2.30|) for i?^ ^ 1 
and <^ 1. First, for the analysis in the region ^ 1 and also for later use, we present 
another expression of t„ obtained by applying the Poisson's resummation formula. 



E 9in)= Yl / dy g{y) e'-^-^y 



to the n^-summations in (j2.3(J|) : 

-1/2 f 1 



= (det g ) exp I 



2(Q-i) 



00 

N-2 ^ N-2 



N-2 

C 



(3.3) 



X E exp < -27r^i?^ ^ m^^Q ^ rUs + 27iin ^ nfm^ > , (3.4) 

mi ,mjv— 2=— oo L r,s=l r=l J 

where the matrix Q is the lower-right (A^ — 2) x (A^ — 2) part of Q, 

Qrs = Qrs, (r,s = l,2,--- ,iV-2), (3.5) 

and is defined by 



)rsQso, (r = l,2,---,iV-2). (3.6) 
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In obtaining the expression (j3.4p . we have used that 



iV-2 



QoO — ^ Qor{Q )rs^'50 



r,s=l 



(3.7) 



00 



which is vahd for any matrix Q and its (A^ — 2) x (A^ — 2) submatrix Q . 

Some comments on the formula ()3.4j) are in order. First, {n^} in ()3.(j|l is nothing but 
the configuration which, without the constraint that be integers and keeping n = 1 fixed, 
minimizes the Hamiltonian (j2.28|) . 



N-2 



H{nr; n) = ^^Q rs^r-ns + Qorrir + ^ Qoon'^ ■ (3.8) 

r,s=l r=l 

Figure 121 shows the configurations {n^} in the cases of 6 = 0.1 (left figure) and 6 = 10 (right 
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-0.1 ■ 
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-0.25 



6 = 10 



Figure 2: The minimum energy non-integer configurations {n^} with n = 1 in the case 
N = 100. The value of 6 is 6 = 0.1 in the left figure and 6 = 10 in the right one. 

figure) for = 100. As seen from the figure, {n^} is localized around r = (mod A^) to 
screen the charge n = 1? 

Our second comment is on the term exp (-[2(Q-i)gJ (^V^^)) in (El- The exponent 
is equal to the value of the Hamiltonian H for the configuration {n ■ n^}, namely, the (non- 
integer) configuration minimizing H for a given n: 

2 

H{n.n';-n)= I" . (3.9) 

As was analyzed in (TH], 1/ (Q~^)oo finite in the limit A^ oo.^ Figure El shows limTv^oo [2(Q~^)qq] ^ 

^As seen from figure |21 the charges near r = all have opposite sign to no for larger values of 6, while 
have alternating signs for smaller h. 

3l/(Q-i)„^ is related to S'^^ in m by l/(Q-')oo = ^ (1/2 + ^^o/(l - ^oo))- 
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Figure 3: liniAr^oo [2((5 ^)qq\ ^ a function Info. The dots represent liniAr^oo [2(Q ^)qq\ ^ 
at 6 = 1/100, 1/10, 1, 10, 100, 1000 obtained by evaluating its values for N = 50, 100, 200, 
300, 400, 500, 600 and then extrapolating them to = cx3 by using the fitting function of 
the form ^l=QCk/N^. The curve interpolating these six points is 0.429214 + 0.216204 x lnfo + 
0.0379729 x (In 6)2 + 0.00186225 x (In 6)3. 

as a function of In b. It is a monotonically increasing function of b and, as we shall see in section 
4, approaches zero as 6 — > 0. 

Now let us consider t„ for i?^ ^ 1. Since Q is a positive definite matrix, the configuration 
rrir = for all r = 1, 2, ■ ■ ■ , — 2 dominates the mr-summation in (j3.4p for ^ 1. Namely, 
we have 

t„~R'^'-^<ietQ) exp - . ^ . (R' » 1) . (3.10) 



This expression can also be obtained by simply replacing the n^-summations in ()2.30p with 
integrations over continuous variables Pr = rir/R. Since the coefficient of in the exponent 
is non-vanishing for 6 > 0, this t„ cannot lead to a desirable rolling profile. 

On the other hand, tn ()2.30p for i?^ ^ 1 can be approximated by 

tn^exp(^-^H{niin);n)'^ , (i?' < 1) , (3.11) 

where {nl.{n)} is the mteg'er configuration which minimizes the Hamiltonian for a given n. This 
configuration {nKn)} is in general different from but is close to {n ■ n^'}, the configuration 
minimizing H{n^] n) without the integer restriction. Therefore, rewriting (j3.1ip as 

t„~exp( --^^^ ) xexpf-^A^nl , (3.12) 



2(Q-i), R^l "V 



'00 



with AH[n) defined by 

AH{n) = H{nl{n); n) - H{n ■ n^; n) , (3.13) 
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the second factor of (j3.12p is expected to have a milder n-dependence than the first one. Figure 
13 shows AH{n) for = 2048 and 6 = 0.1. We see that AH{n) is in fact roughly proportional 
to n. Therefore, we cannot obtain a rolling profile in the case i?^ -C 1 either. 



AH{n) 



50 - 
40 - 
30 - 

20 - 



i 



n 



Figure 4: AH{n) v.s. n for = 2048 and b = 0.1. For obtaining AH{n), we approximate 
the integer- valued charge n^(n) by the integer nearest to n • n^. However, this {nj.(n)} does 
not necessarily satisfy the constraint n + Y^^~^ nj.(n) = 0. In the figure, only the points which 
satisfy the constraint are plotted. Distribution of the points is insensitive to the value of N if 
it is large enough. 



Summarizing this subsection, for both i?^ S> 1 and -R^ <^ 1, the coefficient t„ of the 
component field t{x) has the leading ra-dependence of the form 

Then, t{x) itself shows the behavior 

t{x) ~ exp ^-(Q~^)oo(3;)^^ X (oscillating part) , (3.15) 

and cannot approach the tachyon vacuum as a; ^ oo. This is the case even if we adopt the 
alternating sign solution (— l)"t„. 



3.2 Numerical analysis using Monte Carlo simulation 

For studying t„ and t{x) for intermediate values of R^, we shall carry out Monte Carlo simu- 
lation of the Coulomb system with partition function (j2.27|) . Hamiltonian H and temperature 
R"^. We have adopted the Metropolis algorithm. Since the total charge must be kept zero, a 
new configuration is generated from the old one {n,.} by randomly choosing two points r and 
s on the circle and making the change {nr,ns) —>■ {rir + l,ns — 1). This new configuration is 
accepted/rejected according to the standard Metropolis algorithm. 
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3.2.1 Time derivatives of In t(£c) 



First we investigate the time derivatives of the logarithm of t{x). Defining T{x) by 

t{x) = e^(^) , 

we have 

dT{x) 1 



dx R 
d^T(x) 1 



dx^ i?2 

where the average {O)^ for a given x is defined by 



{{n\ - {n)l) 



-. oo 

(Q) = _L Q -H{nr;n)/R2+nx/R 



n,ni,--- ,njv-i=— oo 



(3.16) 

(3.17) 
(3.18) 



(3.19) 



dT{x)/dx 



X 



Figure 5: The numerical results of dT{x)/dx at x = 0.5, 1.0, 1.5, ■ ■ ■ , 50.0. Here we have taken 
6 = 10, i?^ = 1 and N = 256. These points are well fitted by a linear function 0.4530x — 0.1713, 
and the slope 0.4530 is close to (Q~^)oo = 0.4497 for the present b and A^. The small oscillatory 
behavior can be better observed in d'^T{x)/dx'^ shown in figure El 



i?2 


0.3 


0.5 


1.0 


1.5 


2.0 


3.0 


5.0 


(Q~')oo 


slope 


0.4557 


0.4540 


0.4530 


0.4502 


0.4498 


0.4498 


0.4497 


0.4497 



Table 1: The slope of the linear function of x obtained by fitting the Monte Carlo results of 
{n)x/R with fe = 10 and = 256. They are almost independent of and close to (Q~^)oo- 
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The numerical results of the "velocity" (j3.17|) versus x for b = 10, R"^ = 1 and = 256 
are shown in figure El We find that dT[x)/dx is almost linear in x. The slope of the fitted 
linear function for the various and = 10 and N = 256 as well as the value of (Q~^)oo 
for the present b and are shown in table H These results show that the behavior of ()3.15p 
obtained in the regions R^ <^ 1 and i?^ ^ 1 is valid also in the intermediate region of R^ . 
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Figure 6: (fT{x)/dx'^ at x = 0.5, 1.0, 1.5, ■ ■ ■ ,50.0. Here we have taken 6 = 10, i?^ = 1 and 
= 256. The horizontal line shows the value of (Q~^)oo = 0.4497 for the present b and A^. 



The Monte Carlo results of the "acceleration" (j3.18p at various x are shown in figure IHl 
in the case of 6 = 10, R^ = 1 and A^ = 256. The acceleration oscillates with period roughly 
equal to 8. The center of the oscillation is around {Q~^)oo = 0.4497, which is consistent with 
fj3.15|) .^ We have studied the acceleration for other values of R^ and found that it oscillates 
around (Q^^)oo for any R^. 

3.2.2 Numerical analysis of t„ 

The n-dependence of tn itself can be directly measured using Monte Carlo simulation as follows 
[HHj . Here, we use instead of R^ the inverse temperature /3 = 1/-R^, and make explicit the 
/3-dependence of tn to write it as tn,i3- Let us define the average (C)n,/3 with subscript n and 
/3by 

^ oo 

tn B 

^> ni,--- ,njv-i=— oo 

(« + E^=l'«r = 0) 

^In this case, the parameter a of (|3.1() is a = [2((5~^)oo^^] ^ — 1.112. If the subleading part i„ of H3.1|l is 
independent of n, the oscillation period of t{x) should be 2a — 2.224. The fact that the period of figure El is 
nearly equal 8, which is four times the naive period, suggests that t^k+i, tAk+2 and tik+3, are negligibly small 
compared with t^k- 
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Integrating the relation 
with respect to (3, we obtain 



d 

— liatn,/3 = {-H)n,l3 , 



tn,(i=0 



exp (^j^ d(3'{-H) 



n,f3' , 



(3.21) 



(3.22) 



and hence 



^IhP^exp ( r dp' {{~H)^,f,, - {-H)^=o,p,)] . (3.23) 

„=0,/3=0 \Jo J 

Eq. (j3.1(J|) imphes that t„,/3=o is independent of n and therefore t„^^=o/^n=o,/3=o = 1- Thus, we 
obtain the formula 



in=0,f3 t 



t 



11=0,13 



exp jy(3'{{H)r,,p> - (i^)„=o,/3')) 



(3.24) 



This allows us, in principle, to directly evaluate the n-dependence of t„ using the expectation 
values of H obtained by Monte Carlo simulation. Note that {H)n,/3 in high temperature region 
/3 -C 1 is given using (|3.1(jp by 



(if)„,^ ^ + ■ n% {(3 « 1) . 



2/? 2(g-i; 



(3.25) 
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Figure 7: (-ff)„,/3 - (iJ)„=o,/3 - n^/(2(Q"^)oo) for n = 10 (red points), n = 20 (green), n = 40 
(blue), and n = 100 (purple). Here we have taken N = 256 and b = 10. 



Figure [7| shows the Monte Carlo result of {H)n,f3 — {H)n=o,i3 — '"'^/(2(Q ^)oo)) namely, the 
deviation of {H)n,f3 — {H)n=o,i3 from the high temperature value, for the various values of n 
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{b = 10 and = 256). As (3 is decreased, the data for each n approach zero. On the other 
hand, they seem to approach a positive constant as /3 is increased. Since the asymptotic value 
at large /3 grows no faster than linearly in n, the results of figure [7| together with the formula 
fl3.24|) seems consistent with our low temperature analysis using ()3.12|) and figure El 

From our analyses in this subsection we have found that the behaviors ()3.14|) of tn and 
fl3.15|) of t{x) in the i?^ ^ 1 and i?^ -C 1 regions are valid for any R^. Concerning these 
behaviors, R = 1 does not seem to be a special radius. 



4 Possible rolling solution with 6 = 

Our analysis in the previous section implies that rolling solutions with desirable profile can 
never be obtained unless the coefficient of in the exponent of (|3.14p vanishes. Namely, we 
must have 

= 0. (4.1) 

(g-^)oo 

This condition can in fact be realized by putting 6 = as can be inferred from figure El though 
the precise way how this condition is satisfied differs between the A^ = even and the odd cases. 
For an even (and finite) A^, the condition 1)4.11) is satisfied by simply putting 6 = 0. This is 
because Q with A^ = even and 6 = has a zero-mode {(—1)^}: 

N-2 

^Qrs{-iy = 0, {r = 0,l,--- ,N-2; N = even,b = 0). (4.2) 

s=0 

This zero-mode is at the same time the minimum energy configuration ()3.6|) .^ 

On the other hand, in the A^ = odd case, the condition ()4.1|) is realized by putting 6 = 
and in addition taking the limit N oo. In fact, numerical analysis shows that 

^ ^'^^^+o(^], (AT » 1, AT = odd) . (4.3) 



This property is related to the fact that the matrix Q with 6 = has an approximate zero-mode 
for large and odd A^. This approximate zero-mode, which is at the same time the minimum 
energy configuration ()3.6|) . is given by 

^ (-1)^ X Tl - V (r = 0,±l,±2,---,±(Ar-l)/2), (4.4) 



^Note that the condition for the minimum energy configuration is X^fLo^ Qrsn'^ — for r = 1, 2, • • • , iV — 2, 
while the zero-mode of Q should satisfy this equation for r = 0, 1, • • • , N ~ 2 including r — 0. 
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where the index r should be understood to be defined modiV; This configuration 

satisfies ^1^=0 = and 



A proof of ()4.5j) is given in the appendix. 

In the rest of this section we shall consider only the case of = odd since (j2.12|) is closed 
among the | = odd) states (| A^) is the state given by ()2.20|) in the present case). Here we 
shall just make a comment on the case of A^ = even. When A^ = even and 6 = 0, the coefficient 
tn is independent of n as can be seen from the expression (D with = (-1)^ Therefore, 
we have t{x) = to Yl'^=-ooi^^)'^^'^^^^- ^ naive summation of this series gives t{x) = 0, or the 
Poisson's resummation formula gives t{x) = 2nto J2m=even/odd ^ {i^l R ~ t^t^) EE]- 

We would like to investigate whether the solution with 6 = and A^(= odd) — >■ oo can be 
regarded as a rolling solution. However, it is not an easy matter to repeat the Monte Carlo 
analysis of Sections 13.2.11 and 13.2.21 in the present case. First, since the coefficient (Q~^)oo of 
the leading x-dependent term of ()3.15|) blows up as A^ — >• oo,^ so do the slope of the velocity 
curve of figure El and the central value of the acceleration curve of figure IHl Therefore, it is 
hard to read off the subleading x-dependence which should become the leading one in the 
limit A^ ^ oo. Second, the direct evaluation of the coefficient t„ using the formula (j3.24j) is 
not an easy task because of bad statistics problem. Namely, the difference (-ff)n,/3 — {H)n=o,i3, 
which is of order n^/(2((5"^)oo) (see (j3.25|) ). is very small compared with the leading bulk 
term (A^ - 2)/(2/3) of {H) when 6 = 0. 

Here we shall content ourselves with the analysis of t„ in the low temperature region 
i?^ <^ 1 using the expression ()3.12p . Since the leading term exp(— n^/(2((5"^)oo-R^)) of ()3.12p 
disappears in the limit A^ — oo, we study the n-dependence of AH{n) in the same way as we 
did in figure El The result for A^ = 8191 is plotted in figure |H1 

Note that the data of figure IHlhas an approximate periodic structure with respect to n with 
period of about 4100. This periodicity can be understood from ()3.4j) for tn and the expression 
()4.4j) for n^' independently of the low temperature approximation. In fact, ()3.4|) without the 
leading term exp(— ra^/ (2(Q~^)oo-R^)) is invariant under the shift n ^ n + {N + l)/2 since the 
change of the exponent in the m^-summations under this shift is an integer multiple of 27ri for 

of (j4.4j) . Although figure |H1 shows data only for positive values of n, recall that t„ is even 
under n —n; (j2.29|) . This parity symmetry and the periodicity lead to the structure shown 
in figure |H1 

This divergence is only an apparent one coming from applying the evaluation of t{x) given in to the 
case with an infinitesimally small a. If the leading term H3.14|l of tn is missing from the start, we have to adopt 
a different way of estimating the summation (|2.26() . 




(4.5) 
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Figure 8: AH{n) v.s. n for iV = 8191 and 6 = 0. 



The periodicity stated above, 

tn+{N+l)/2 = tn , (4.6) 

is not an exact one for a finite since both the condition 1)4.11) and eq. ()4.4p for n'^ are 
only approximately satisfied. If the periodicity ()4.(ip were exact, t{x) given by ()2.2(ij) could be 
rewritten as 

t (a;) = tone-period (a;) ^ exp ( kx\ , (4.7) 

fc=-oo ^ ^ 

where tone-period (a;) is defined by 



[(7V+l)/4] 

tone-period (a^) ^ ^ t^ 6 ^ , (4.8) 

n=-[N/A] 

with [c] being the largest integer not exceeding c. Namely, tone-period (a;) is the one period part 
around n = in the summation of t{x). The geometric series multiplying ()4.7|) is formally 
summed up to give an unwelcome result; it is equal to zero or the summation of delta functions 
for pure imaginary values of x. However, since the periodicity ()4.6j) is not exact for a finite A^, 
we have to carry out more precise analysis taking into account the violation of the periodicity 
to obtain the profile t{x) in the limit oo. Here we would like to propose another way of 

defining t(x) which could lead to a desirable rolling profile. It is the N ^ oo limit of the one 
period summation ()4.8p : 

t(x) = lim tone-period (a;) • (4.9) 

This is also formally equal to the original summation (|2.26|) with N = oo. 
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Let us return to the low temperature analysis with i?^ ^ 1. For t{x) given by (j4.9p . it is 
sufficient to study the ra-dependence of AH{n) only in the half period region n G [O, [(A^ + 
l)/4]] , which is shown in figurelHlfor = 8191. As shown in figure|3 AH{n) has a complicated 




Figure 9: AH{n) for n G [0,2048] in the case A^ = 8191. The red line is an auxiliary one 
with slope 3/2. 



structure with peaks and valleys. However, if we neglect such local structures and see figure 
ini globally, we find that AH{n) grows almost linearly in n; AH{n) oc n (the red line in figure 
El is the line with slope 3/2). This could give a desirable t„ with the behavior ()1.4j) up to 
complicated local structures. However, it is a nontrivial problem whether the N ^ oo limit of 
(14. 9j) really exists even when we take into account the local structures. Here we shall point out 
a kind of self-similarity of AH{n) and hence of t„; 2 AH{n)\]y ^ AH{2n)\2N- Figure ITUl shows 
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Figure 10: AH{n) for A^ = 8191 (blue points) and that for A^ = 4095 (red points). Both the 
horizontal and the vertical scales are doubled for the red points. 
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AH{n) for = 8191 (blue points) and that for N = 4095 (red points). Both the horizontal 
and the vertical scales are doubled for the N = 4095 points. For example, the real coordinate 
of the red point (1500, 1600) in the figure is actually (750, 800). Note that the red points and 
the blue ones have overlapping local structures. It is our future subject to study whether t{x) 
of ()4.9p can exist for tn with such self-similarity. 



5 Thermodynamic properties of t(x) 

In the BCFT analysis, the rolling tachyon solution is obtained by the inverse Wick rotation 
of one space direction compactified on a circle at the self dual radius R = 1 [2HIC]- Therefore, 
also in our VSFT construction of the rolling tachyon solution, it is natural to expect that the 
meaningful solution exists only at R = 1. This is also supported by the fact that the correct 
tachyon mass squared —1 has been successfully reproduced in the analysis of the fluctuation 
modes around the D25-brane solution in VSFT j23 12^1 If the tachyon mass squared is 
equal to —1, the natural mode of the expansion in ()2.26|) is e"'^ with R = 1 [31111; in particular, 
the n = ±1 modes e^^ are the massless modes at the unstable vacuum. 

In this section we shall study how this critical radius R = 1 appears in our construction of 
time-dependent solution in VSFT, especially in the component field t{x). One would naively 
expect that the N oo limit of our solution ()2.24|) can exist only at i? = 1. Here, we do not 
pursue this possibility directly, but instead address the problem from a statistical mechanics 
point of view. Recall that t{x) (j2.26p and t„ (j2.27|) can be interpreted as the partition functions 
of a statistical system of charges located on a unit circle with temperature R^. One possible 
mechanism of i?^ = 1 being a special point for this statistical system is that it undergoes some 
kind of phase transition at i?^ = 1. In this section, we first claim, on the basis of a simple 
energy-entropy argument, the presence of a phase transition at i?^ = 1. Then, we study in 
more detail the thermodynamic properties of the system both analytically and numerically. 
Our results here suggest but do not definitely confirm the presence a phase transition at 
R^ = 1. 



5.1 Boundstate phase and dissociated state phase 

Let us consider the statistical system with partition function t{x = 0)-7 

m= exp(-^H{nr;no)Y (5.1) 

no,ni,--- ,njv_i=-oo ^ ^ 

(no+niH hnjv-l=0) 

^Although we consider here t{0) with x = 0, t{x ^ 0) and <„ have the same bulk thermodynamic properties 
since the difference is only the local one at r = 0. 
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We would like to argue that this system has a possible phase transition at temperature i?^ = 1. 
In the low temperature region i?^ -C 1, configurations with lower energy contribute more to the 
partition function. The lowest energy configuration of the Hamiltonian H ()2.28|) is of course 
that with all = 0. Due to the self-energy part 2 ln(A^/4) of Qrr ()2.25|) . the energy of a generic 
configuration with zero total charge can be In A^-divergent. Finite energy configurations are 
those where the charges are confined in finite size regions and the sum of charges in each 
region is equal to zero. Namely, they consist of neutral boundstates of charges. The simplest 
among them is the configuration of a pair of +1 and —1 charges with a finite separation. Let 
us consider the configuration {nl^''^'*} with = +1, n^+A = — 1 and all other = for a 
given position k and separation A. The energy of this configuration is 



21n( ^ 1 + 21n 



2 sin 



TT 



A 



N 



(5.2) 



This is approximated in the close case A 
A = O(A^) (mod N) by 



0(1) (mod A^) and in the far separated case 



- + 2 In — ^ 
2 V 2 

21nA^, 



A 
A 



0(1), 
0(iV). 



(5.3) 



The energy in the finite separation case A = 0(1) is indeed free from the InA^ divergence. 
Other neutral boundstates with finite energy are, for example, a pair of charges +q and —q 
with g > 1, and a chain of alternating charges 



(0,0- ■• ,0,q,-q,q,-q,--- ,q,-q,0, 



.0). 



(5.4) 



Configurations with isolated charges have In A^-divergent energy as seen from the A = 
0{N) case of (j5.3j) . However, this does not imply that such configurations do not contribute 
at all to the partition function: we have to take into account their entropy. Let us consider 
a naive free energy argument of an isolated charge. As seen from Qrs (j2.25|) or (j5.3|) for 
A = 0{N), the energy of an isolated charge is £^ = InA^, while the entropy of this charge 
is S = In A^ since there are A^ points where it can sit on. Therefore the free energy of this 
isolated charge is given by 

^ isolated charge = ^ — = (1 ~ ^ ■ (5-5) 

This means that, for i?^ > 1, the free energy becomes lower as more isolated charges are 
excited. Namely, R"^ = 1 could be a phase transition point separating the boundstate phase 
in i?^ < 1 and the dissociated state phase in i?^ > 1. To confirm the existence of this phase 
transition, more precise analysis is of course necessary. 
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5.2 Dilute pair approximation 



Before carrying out numerical studies of the system (|5.ip for the possible phase transition at 
i?^ = 1, we shall in this subsection present some analytic results valid in low temperature. In 
the low temperature region i?^ -C 1, it should be a good approximation to take into account 
only the pairs of charges as configurations contributing to the partition function ()5.1|) . This 
approximation is better for larger b since more complicated boundstates such as ()5.4j) have 
larger energy coming from the 6/2 term of Qrr ()2.25|) . The partition function of a pair of 
charges is given by summing over the position k and the separation A of the configuration 

N-l N-1 

-H({n(^--^'})/R2 



^1- 



pair 



N-1 N-1 
k=0 A=l 



(5.6) 



In the low temperature region where the number of pair excitations is small and the pairs are 
far separated from each other, we can exponentiate 2'i_pair to obtain 

m = exp (Zi.pair) . (5.7) 

We call this "dilute pair approximation". Using ()5.2|) . Zi.p^h is calculated as follows: 

-2/_R2 



pair 



iV-1 N-1 

EE' 

k=0 A=l 



-b/{2B?) 





-2/_R2 


ttA 












AT 



-2/_R2 



2N 



2/i?2 _ 1 



vre 



6/4 



^2 rn{l-l/N) 

J-n/N 
-2/_R2 



-2/_R2 



dy{smy) "^^^^ 
{R^ < 2) 
(i?2 > 2) 



(5.8) 



2; v^r(i-^) 

In the second line of fl5.8j) we changed the A-summation to the integral with respect to 
y = ttA/N for ^ 1. This integral is divergent (convergent) at the edges as — > 00 in the 
region R^ < 2 [R'^ > 2). The final result of (j5.8p in i?^ < 2 has been obtained by taking the 
contribution near the edges y rr/N and 7r(l — 1/A^), while that in i?^ > 2 by extending the 
^-integration region to [0, vr]. 

As we raise R^, the more pairs are excited and the separation of the two charges of a pair 
becomes larger, leading to the breakdown of the dilute pair approximation. From ()5.8p . one 
might think that this breakdown occurs at i?^ = 2. However, let us estimate more precisely 
the limiting temperature above which the dilute pair approximation is no longer valid. The 
criterion for the validity of the dilute pair approximation is given by 



{p)xA<N., 



(5.9) 
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where (p) and A are the average number of pairs and the average separation between +1 and 
the —1 charges of a pair, respectively. First, (p) is given by 



The average separation A is calculated as follows: 



pair 



(5.10) 



IT N 

A = — X 

TT 



7r(l-l/Ar) 



dy min(?/, vr — y) (siny) '^^^ 



tt/N 



tt{1-1/N) 



dy {sin y) '^/^ 



tt/N 



1 - i?V2 



(i?2 < 1) 



1 



v^r(i-^) 



(5.11) 



(1< i?2 < 2) 



l-i?2 ' 

' ' 'r(i-^) 

where the separation is the smaller one between A and — A. We are considering only the 
region < 2 since we have (p) ~ N^^'^/^'' ^ in the other region R^ > 2. Note that the 
critical R^ below which the y-integration in the numerator of ()5.11|) diverges at the edges has 
changed to i?^ = 1 due to the presence of mm{y, n — y). 

As seen from fl5.1ip . the average separation A diverges as R^ | 1. This is consistent with 
the boundstate-dissociated-state transition which we claimed to occur at i?^ = 1. From ()5.10|1 . 
(|5.8|) and (|5.11|) . the condition (|5.9|) for the validity of the dilute pair approximation (j5.7j) is 
now explicitly given by 



1/R^ 



Tce 



6/4 



<1, (^'<1) 



(5.12) 



The breakdown temperature of the dilute pair approximation determined by ()5.12|1 becomes 
lower as b is decreased. Moreover, for smaller b we have to take into account also longer neutral 
boundstates such as ()5.4jl . and hence the dilute pair approximation becomes worse than the 
above estimate. However, we would like to emphasize that the free energy argument of an 
isolated charge using ()5.5|) holds independently of the value of b. 



5.3 Monte Carlo analysis of internal energy and specific heat 

In order to study whether the boundstate-dissociated-state phase transition which we predicted 
in section 15.11 really exists, we have calculated using Monte Carlo method the internal energy 
E and the specific heat Cy of the system t{x = 0): 

^ = -|j^ln«(0). (5.13) 
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(9^ 1 



(5.14) 



Figures ITT] and IT^ show E and Cy for b = 10 and & = 1, respectively (A^ = 512 in both 



E 

2 - 




-4 — » — » — t — » 



Figure 11: E and Cy vs. i?^ for N = 512, 6 = 10 and x = 0. The curves in the smaller 
region of and the straight lines in the larger region have been obtained by the dilute 
pair approximation and the high temperature approximation, respectively. 
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Figure 12: E and Cy vs. R^ for = 512, 6 = 1, and x = 0. 

the figures). The curves in the smaller R^ region have been obtained by using the dilute pair 
approximation ()5.7|) : 



2 /vre^/^^ 



— lnt(0)~ ^ 
N ^ ' 2(3-1 \ 2 



{R^ < 1) 



(5.15) 



The staright lines in the larger R^ region are from the high temperature approximation (c.f. 
(jSini)): 

i_lnt(0)^-iln/5, (i?2»l), (5.16) 
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giving 

E^\r\ Cv^\. (5.17) 

The curves of the dilute pair approximation fit better with the data in the 6 = 10 case than 
in the 6 = 1 one. This is consistent with our analysis in section EIH using ()5.12|] . 

Figures ^2 and ^1 show no sign of phase transition around = 1. Note that there is a 
peak strucrure in Cy- The of the peak becomes larger as the parameter h is increased. We 
have carried out simulations for smaller values of A^, and found that Cy, and, in particular, 
the height of the peak are almost indepdendent of A^. Therefore, the peak in Cy is not a 
signal of a second order phase transition.^ 



E Cv 




Figure 13: E and Cy for 6 = and = 511. Here, the horizontal axis is (3. The curve of 
E in the dilute pair approximation in the larger f3 region is invisible since it almost overlaps 
with the /5-axis. 

We claimed in section 0] that we have to set 6 = and = odd for obtaining a sensible 
rolling solution in our construction. Figure IT^ shows E and Cy of t{0) for 6 = and A^ = 511. 
The high temperature approximation ()5.17p fits well with the data in the region ^ 3. 
However, the dilute pair approximation is not a good approximation even in the region of 
the largest f3 in the figure. As we mentioned in section 15. 2^ we have to take into account 
longer neutral boundstates besides the simple pair for obtaining a better low temperature 
approximation.^ In any case, we cannot observe any sign of lower order phase transitions from 
the figure. 

®The global structure of Cy with a peak and the asymptotic value of 1/2 can roughly be reproduced from a 
simple partititon function ^'^^_^exp (— /36n^/4) neglecting the Coulomb interactions although the position 
and the height of the peak differ from those obtained by simulations. 

^In fact, incorporation of the chain excitation H5.4|l with q — 1, whose energy is given by H5.19|l in the next 
subsection for sufficiently large length A, seems to considerably improve the low temperature approximation 
for 5 = 0. 
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5.4 Correlation function 



Next we shall investigate the correlation functions of rir in the low and the high tempareture 
regions. Here, we consider the two-point correlation function {rir Ur+A) in the system with 
partition function t[x = 0) for a large distance |A| ^ 1 (mod A^): 

^ oo 

(rikrik+A) = j7^ ^fc^fc+Ae-^^^"-"") . (5.18) 

^ ^ no,--- ,"jv-i=— oo 
(noH hrtjv-i=0) 

Let us consider the low temperature region i?^ -C 1. There are two candidate configurations 
with lower energy which mainly contribute to ()5.18|1 . One is the configuration of a pair 
excitation, {±ni^''^^}, which we defined in section I^TT] (rik^nk^^) = (±1, =f1) and all the other 
Ur = 0. The other configuration is the chain of alternating sign charges ()5.4|) with ends at 
r = k and k + A; (n^, Uk+i, ■ ■ ■ , Uk+A) = ±(1, —1, ■ ■ ■ , —1, 1) and all other rir = 0. We call 
it the chain excitation. This chain excitation exists only in the case of odd A. In the case 
of even A, there are similar alternating sign configurations with zero total charge. Here, for 
simplicity, we consider only the case of odd A. 

The energy of the pair excitation is given by ()5.2|1 . We calculated numerically the energy of 
the chain excitation i^chain- The dots in figure ITU show i^chain for the various lengths A (A^ = 
1023 and 6 = 0). These points are well fitted by the curve i/chain = 3.42 + 0.500 ln[sin(7rA/A^)]. 




Figure 14: The energy of the chain excitation at various A in the case of = 1023 and 6 = 
(dots). The curve is i^chain = 3.42 + 0.500 ln[sin(7rA/A^)] obtained by fitting. 

We have carried out this analysis for various values of A^, and found that the dependence of 
-f^chain on (sufficiently large) A and A^ is given by 

^chainlfe^o = ^lnN + ^\n 

In particular, the coefficient of lnsin(7rA/A^) is 1/2 and independent of A^. In the case of 
6 7^ 0, Hchain has an additional self-energy contribution (6/4)(A + 1). Thus we find that the 



sm 



TT 



A 



N 



+ const. 



(5.19) 
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contributions of the two kinds of excitations to the correlation function are 



sm 



vrA 



^-l/(2fl2)g-6A/(4ij2) 



sm 



vrA 



-1/(2R2) 



: pair excitation , (5.20) 
: chain excitation . (5.21) 



Eqs. (I5.2(jp and (j5.2H) are vahd in a wide range of A including both A = 0{N/2) and 1 ^ 
A ^ N/2. In particular, in the region 1 ^ A ^ N/2, the A^-dependence cancels out and we 
obtain 



1 



{rik Uk+A) 



l<A<CAf/2 



A2/«2 
-bA/{4R2)_ 



: pair excitation , 
: chain excitation . 



(5.22) 
(5.23) 



AV(2fi2) 

From ()5.20|) - ()5.23p we see the foUowings. In the case of 6 7^ 0, (n^rifc+A) is given by the con- 
tribution from the pair excitation, ()5.20|1 and ()5.22|1 . for a large distance A ^ 1. Contribution 
of the chain excitation is suppressed by e~^'^/*^^^^\ On the other hand, in the case of 6 = 0, 
(ukUk+A) is given by the contribution from the chain excitation, ()5.21|1 and ()5.23|) . 

Next, let us consider the high temperature region ^ 1. For ^ 1, approximating the 
summations over in ()5.18p by integrations over continuous variables Pr = Ur/R, we obtain 



{rikrik+A) -^R\Q 
Incidentally, (Q^^)^ with A; = is related to ()3.6p by^° 



(5.24) 
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detQ 
detQ 



n 



c 



(5.25) 



We calculated {Q~^)k,k+A numerically and found that it is independent of the position k; 
namely, the translational invariance holds for large N despite that the self-energy 6/4 is missing 
for uq (see (I2.25j) ). Therefore, (Q"^)^^.,.^ is essentially equal to up to a A-independent 
factor. 

The A-dependence of is shown in figure |21 in the cases of b = 0.1 and 10. As we 
mentioned in the footnote there, and hence {Q~^)k,k+A has quite different A-dependences 
in the smaller A region: for larger b, {Q~^)k,k+A is negative definite, while it has alternating 
sign structure for smaller b. However, for A = 0{N/2) in the mid region, {Q~^)k,k+A is 
negative definite and has the following universal A-dependence for any non- vanishing b as we 
shall see below: 



{Q 



-1^ 



)k,k+A 



N 



-2 



sm 
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(A = 0(A^/2)) 



(5.26) 
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Figure 15: ln|(Q ^)oa| for = 511 and 6 = 10 (dots) and the fitted red curve —12.02 — 
2.1071n|sin(7rA/A^)|. 

First, for larger b, {Q^^)k,k+A is well fitted by (j5.26|) for any A (see figure HHI). We have 
confirmed the A^-dependence of (I5.26p by the fitting for various A^. Next, let us consider 
{Q~^)k,k+A for a smaller b. Figure ITBl shows {Q~^)oa in the case of 6 = 0.01 and = 511 
in all the region of A (left figure), and In |(Q~^)oa| in the mid region 180 < A < 332 where 
(Q~^)p^ is negative definite (right figure). As we see from the right figure, the A-dependence 
of ()5.26|) holds well in the mid region of A. By carrying out this analysis for various A^, we 
have confirmed the A^-dependence of fj5.26|) also for smaller b. 
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Figure 16: The left figure shows (Q ^)q^ for b = 0.01 and A^ = 511 in all the regi 



ion 



A = 1, - ■■ ,A^ — 1. The right figure shows In | (Q ^)q^I ^-nd the fitted red curve —13.09 — 
2.001n|sin(7rA/A^)| only in the mid region 180 < A < 332. 



The behavior of (Q ^)q^ in the case of 6 = and A^ = odd is quite different from the 
non-zero b case above, in this case is given by ()4.4p . and {Q~^)q^ is shown in figure IT7I for 



-"^^ Numerical analysis shows that det Q / det Q is finite in the limit N 
det g / det Q ~ 0.59 TV for 6 = 0. 



oo for 6 7^ 0, while we have 
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= 511 and 6 = 0. {Q ^)oa is a linear function of A with alternating sign for all A. 
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Figure 17: {Q'^)q^ for N = 511 and 6 = 0. 



From the above analysis of {ukUk+A) in both the <^ 1 and i?^ ^ 1 regions, we find the 
foUowings. In the case of 6 7^ 0, {rik n^+A) is given by (j5.20|) in the i?^ <^ 1 region, and by (j5.26|) 
m the > 1 

region. Namely, {n^nk+A) has the same kind of A-dependence | sin(7rA/A^)| 
although the exponent 7 differs in the two regions. On the other hand, 't-zc+a) in the 
6 = case has entirely different A-dependences in the two regions; it is given by (j5.2H) in 
the <^ 1 region, and by (j4.4p in the ^ 1 region. This result suggests the existence 
of a phase transition at an intermediate R^ at least in the 6 = case. Further study of the 
correlation function in the intermediate region of R^ using Monte Carlo simulation is needed. 

6 Summary and Discussions 

In this paper, we constructed a time-dependent solution in VSFT and studied whether it can 
represent the rolling tachyon process. Our solution is given as the inverse Wick rotation of 
the lump solution on a circle of raduis R which is given as an infinite number of ^-products 
Qb*Qf,* ■ ■ ■ *Qh of a. string field Qb with Gaussian momentum dependence e"''^^/^. We focused 
on one particular component field in the solution, which has an interpretation as the partition 
function of a Coulomb system on a circle with temperature R^. Our finding in this paper is 
that, for the solution not to diverge in the large time limit, we have to put 6 = and take 
the number of Qb consituting the solution to infinity by keeping it even. We also examined 
the various thermodynamic quantities of our solution as a Coulomb system to see whether the 
self-dual radius R = 1 has a special meaning. We pointed out a possibility that R = 1 is a 
phase transition point separating the boundstate phase and the dissociated state phase. Our 
analysis of the correlation function for 6 = supports the existence of the phase transition. 
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Many parts of this paper are still premature and need further study. The most important 
among them is to study in more detail our solution with 6 = 0: whether the limit iV — »• oo 
really exists, and if so, what the profile will be. In this paper we tried to find a special nature 
of the critical radius i? = 1 in the thermodynamic properties of the system. However, we 
have to find a more direct relevance of = 1 to our solution. For example, the most natural 
scenario is that the limit N —>■ oo can exist only at i? = 1. Analysis of the whole of our 
solution not restricted to the component t{x^) is also necessary. 

Originally our time-dependent solution had two parameters, b and R. If we have to put 
6 = and R = 1 for obtaining a solution with a desirable rolling profile, this solution seems 
to have no free parameters. However, the rolling solution should have one free parameter 
which corresponds to the initial tachyon value at x° = 0. It is our another problem to find the 
origin of this parameter. It might be necessary to generalize our solution to incorporate this 
parameter. 

After establising the rolling solution, our next task is of course to apply our solution to 
the analysis of unresolved problems in the rolling tachyon physics. 
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A Proof of eq. (14.51) 



In this appendix we present a proof of (|4.5|) for given by the RHS of (j4.4j) . Since this 
satisfies = 0, and for such we have 

Af-2 N-1 N-1 

XI '^rs^ = XI ^'-^^^ ~ '^JV-l.s^ , (A.l) 

s=0 s=0 s=0 

eq. (j4.5p holds if we can show that 

N-1 . ^ . 

X Qrs^f = (r-independent term) , (r = 0, 1, ■ ■ ■ , iV - 1) . (A.2) 

s=0 ^ ^ 

Before starting the proof of ()A.2j) for of ()4.4|1 . we shall mention eq. ()4.2|1 in the 6 = 
and = even case. This fl4.2p holds owing to a stronger equation ^fljo^ Qrs{—^Y = 0, which 
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is rewritten explicitly as 

N-l 



i;,-i)M„|2.„(^)|^..(^) 

r=l ^ ^ 



(A^ = even) . 



(A.3) 



Formulas essentially equivalent to ()A.3|) can be found in the various tables of series and 
products. 

Now let us consider the LHS of ()A.2|) for given by the RHS of (j4.4j) . Taylor expanding 
ln|2sin[7r(r - s)/A^]| = Reln(l - e^'^*^^-^)/^) in Qrs (IT^ in power series of e^'^*^^-^)/^, we 
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(A.4) 



with /fc defined by 

'^'^ ~ N+1 (^^iTTk/N _^ g-i7rfc/Ar^2 

In obtaining the last line of ()A.4|) we have applied the formula 

M-l 



^ ^ ' M J m' 



s=-A/+l 



1-^ 



iV+l 



to the case of M = (A^ + l)/2 and z = -e^^^^/^, and used that (-£±2'^*'=/^) 
Note that fk defined above has the periodicity: 

fk+N = fk ■ 



(A.6) 



-1) 2 X 



(A.7) 



Introducing the cutoff LN in the fc-summation in ()A.4|) and using the periodicity to make the 
manipulation Y.k=i fk/k = XlfcLi fk Z]p=o 1/ (P^ + k),we have 
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L>1 



1 ^ 



k=l 



InL-ibi — 



2 |r| 

ivTI 



[In(LAr) + 7] 



(A.8) 



where ijj{z) is the polygamma function: 

d 
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dz 



^^\n + z n + lj 



(A.9) 



In obtaining the last hne of ()A.8|) . we have used ipi'^) = ~1 and the asymptotic behavior of 
ipiz) for \z\ ^ 1: 



tp{z) ~ \n{z — 1) + 



2{z-l) 12(^-1)2 
Now we have to carry out the two summations in ()A.8|1 . 
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for large A^. One way to evaluate S\ is to approximate it by a contour integration with respect 

1 r fjz 1 2\r\ 

Si=^— i -/(^) = Res-/(z) = l--^, (A.ll) 
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where f{z) is 



-z (-z)^ + (-l/^)-f--2 
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(A.12) 



A^ + 1 (l + zY 

Note that f{z)/z is regular at z = —1. 

Another way of evaluating Si is to observe the foUowings. fk is of 0{1/N) except at 
k ~ (A^ + l)/2 where fk = 0{N). Therefore, we have only to carry out the fc-summation 
around k ~ (A^+l)/2. Expressing k ask = (A^+l)/2+£, fk is expanded around k = {N+l)/2 
as 



fk={N+i)/2+e 
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and we regain the same result as ()A.11|) : 
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where we have used the formulas, 
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The other summation 5*2 is evaluated in the same manner: 
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where we have used that = —7 — In 4, and that the summations 
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vanish due to the oddness under — 1. Recalling that 
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(A.15) 
(A.16) 



(A.17) 



(A.18) 



and plugging the results (IA.14j) (or (lA-lljl ) for 5*1 and (IA.17jl for S2 into (lA.lSj) . we find that 
the terms on the RHS just cancel out. Therefore, the r-independent constant on the RHS of 
is in fact equal to zero. 
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